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Abstract 

The nature of the phase transition in the lattice Gross-Neveu model with Wilson fermions 
is investigated using a new analytical technique. This involves a new type of weak coupling 
expansion which focuses on the partition function zeroes of the model. Its application to the 
single flavour Gross-Neveu model yields a phase diagram whose structure is consistent with that 
predicted from a saddle point approach. The existence of an Aoki phase is confirmed and its 
width in the weakly coupled region is determined. Parity, rather than chiral symmetry breaking 
naturally emerges as the driving mechanism for the phase transition. 



1 Introduction 



In continuum QCD, the conventional explanation for the smallness of the mass of pseudoscalar 
vr mesons is the following: QCD with Nf massless quark flavours has a global chiral U{Nf) x 
U{Nf) symmetry, which, spontaneously broken, reduces to U{Nf) and yields A^^^ Goldstone bosons. 
Explicit breaking of the original chiral symmetry by a small quark mass renders these Nj Goldstone 
bosons massive with correspondingly small mass. To accord with nature, one of these Goldstone 
bosons (the ?7-particle in the case Nj = 2) has to acquire an an additional mass. To explain this was 
known as the U{1) problem in the continuum. Its resolution there comes from the axial anomaly, 
whereby the axial symmetry corresponding to the U{1) subgroup of U{Nf) is explicitly broken by 
a quantum effect, reducing the number of Goldstone bosons to Nj — 1. 

Naive lattice regularization of such a fermionic theory is hindered by the doubling problem, namely 
that a return to the continuum manifests too many fermionic degrees of freedom. This doubling 
problem is resolved by the usage of Wilson fermions. However the extra Wilson term that removes 
the fermion doublers breaks chiral symmetry explicitly. This effect can be traced back to the 
existence of the axial anomaly in the continuum. For this reason the staggered format has often 
been the favoured one for the study of models with chiral symmetry breaking . 

The Wilson action for free fermions in terms of dimensionless fermionic fields -0(71,) defined at the 
sites n of a d-dimensional lattice is 

sP[ip,il^] = - \ ^[i'in)ir - j^,)ipin + + i'in + n){r + -ff,)'ilj{n)] , (1.1) 

where 

1/2k = MQ + dr . (1.2) 

Here, k is the hopping parameter, r is the Wilson parameter, a is the lattice spacing and Mq is a 
dimensionless fermion bare mass parameter. We use d = 2 and r = 1 throughout. This free fermion 
model is the weak coupling limit of an interactive theory in which a bare parameter g measures the 
coupling of the free theory to some interaction. Even if Mq = 0, now, the Wilson term contributes 
to the hopping parameter and there is no obvious chiral symmetry. The question arises - what is 
the status of the chiral phase transition and the ^7(1) problem on the lattice? 

Despite the lack of an obvious chiral symmetry, there exists a host of numerical and analytical 
evidence for the existence of massless pions in the lattice formulation of QCD. These are believed 
to exist on a critical line Kc{g). In the literature, there are two explanations for the existence of 
the critical line and the masslessness of the lattice pions. 

The first of these was sometimes referred to as the conventional explanation Although there 
is no obvious chiral symmetry at non-zero r, the conventional explanation suggests that tuning k 
effects its recovery in some unknown way. Now, with chiral symmetry recovered at Kc{g), the same 
arguments as in the continuum may be applied. 

The second explanation was first forwarded in 1984 by Aoki j^. Here it is accepted that since 
there is no chiral symmetry in the lattice formulation of QCD, its spontaneous breaking cannot 



be responsible for the masslessness of pions. Instead there is an Ising-hke second order parity 
breaking phase transition. In the single flavour case the order parameter for parity symmetry is 
ifji'j^ijj, the operator corresponding to the single vr meson. The parity symmetric phase is where 
(vr) = (ipi^^'ip) = 0. there is also a phase with long range order where (vr) ^ 0. At the transition 
between these phases, a correlation length ^ diverges. This correlation length is identified as the 
inverse of the pion mass, which, hence, becomes zero on the phase boundary. Thus the pion is 
not a Goldstone boson in the Wilson lattice formulation. Aoki also recovered the current algebra 
relation between pion and quark mass (m^ ~ ruq ~ (k — Kc)) by considering the effective meson 
theory as a scalar field theory in four dimensions with mean field like critical behaviour. In the 
multiflavour case the parity symmetry breaking is accompanied by a flavour symmetry breaking 
and, with it, Goldstone bosons in the form of the charged pions. The rj remains massive according 
to Aoki's analysis, and the U{1) problem in the lattice successfully resolved ||3|, ^. 

Two main features distinguish Aoki's QCD phase diagram from the conventional one. Firstly, the 
existence of the phase transition in Aoki's scenario is due to parity symmetry breaking as opposed 
to chiral symmetry breaking in the conventional picture. The order parameter is (ipij^ip) rather 
than (V'V') 01- Secondly, instead of a single critical line extending from the strongly coupled limit 
g = cxD to K = l/2d in the weakly coupled limit g = 0, Aoki's picture involves the existence of two 
such lines and (in QCD) flve critical points linked by four cusps in the weakly coupled zone. 

Aoki's QCD phase diagram is based on inflnite-volume analyses in the limits of strong and weak 
coupling and on an analogy to the Gross-Neveu model, which, except for confinement, has features 
similar to QCD. One of these features is asymptotic freedom, so that in the Gross-Neveu model, 
as in QCD, the continuum limit is taken in the weakly coupled regime. Aoki's scenario in the 
Gross-Neveu model again involves two critical lines spanning the full coupling range, with three 
critical points at zero coupling, linked by two cusps. This picture is based on saddle point methods 
[^. There exists substantial evidence in support of this scenario in the strongly coupled regime 
[P, 1^, ^, 0, ^, 0, 1^. In the weakly coupled regime, however, the evidence has been 

clear cut Q, |5| and this is the region where our attention is focused. Recently, also, Creutz [ p^ ] 
has posed a question as to the size of the Aoki phase. This question is whether the Aoki phase is 
"squeezed out" between the arms of the cusps at non-zero coupling or whether it only vanishes in 
the weak coupling limit 0, |l^ . 

This sets the twofold motivation for this paper. Firstly, a new type of weak coupling expansion is 
developed ||T^. From it, the partition function zeroes of Wilson fermionic models can be extracted 
in a natural way. This weak coupling technique is then applied to the Gross-Neveu model, where 
the existence of an Aoki phase was first suggested . We confine our attention to the single fiavour 
Gross-Neveu model and variants thereof. We also address the question of the "squeezing out" of the 
Aoki phase at weak coupling. This multiplicative approach to the single flavour Gross-Neveu model, 
shows that the width of the central Aoki cusp is 0{g^) while the Aoki phase has not yet emerged 
at this order from the left and right extremes. Furthermore, that parity symmetry breaking is the 
phase transition mechanism emerges in a very transparent way. 



2 The Gross-Neveu Model 



The original motivation for the introduction of the Gross-Neveu model in the continuum [^] was to 
study a renormalisable quantum field theory involving dynamical spontaneous symmetry breaking. 
Such models evolved from four dimensional four-fermi models studied by Nambu and Jona-Lasinio 



[18 1 and are essentially their two dimensional equivalents. The Gross-Neveu model is, however, 
renormalisable and asymptotically free. It is a model of fermions only, which interact through a 
short range quartic interaction. We start with a generalized Gross-Neveu model, whose action, in 
euclidean continuum space, is given by 

= / + M) ^(x) - ^ {i>{x)^{x)f - f (V;(x)i75V'(^))'| , (2.1) 

where 75 = i~^^i^2 and ip{x) is a 2 component fermion field. Note that we have allowed for two 
different four- fermion couplings. This allows for some flexibility to tune in or out the continuous 
chiral symmetry present in the continuum action [|l^, pC| ]. 

We use the following representation for the Dirac 7-matrices in two dimensions, 
so that the chirality operator is 



75 = i Si72 = ( 1 • (2-3) 




Each term in the action (^]^) is invariant under the continuous global U{1) symmetry 

exp {ia)ip{x) , 'ip{x) exp {—ia)-ip{x) . (2-4) 



If, further, the fermion mass M vanishes, the action (2.1) is also invariant under a discrete global 
chiral Z2 transformation 

V'(x) 'Jsi^ix) , Tpix) V'(x)75 . (2.5) 

This is the symmetry of the original (standard) version of the model, in which the last term of 



(2.1) is absent (i.e., = 0). Finally, if the four fermi couplings are tuned such that = Qt^, the 



discrete chiral symmetry is promoted to a continuous one, namely, 

V'(x) — > exp (i075)^(x) , V'(x) ^ V'(x) exp (26*75) . (2-6) 

This cannot be broken since there are no Goldstone bosons in two dimensions due to the Mermin- 



Wagner theorem [21|. Nonetheless, a topological long range order of the Kosterlitz-Thouless type 
could exist in the model |^2|. The Mermin- Wagner theorem refers only to continuous symmetries 
and does not preclude the spontaneous breaking of a discrete symmetry in two dimensions. In 
the continuum Gross-Neveu model, the spontaneous breaking of the discrete 75 symmetry leads 



to dynamical fermion mass generation. The mass term explicitly breaks chiral symmetry and is 
analogous to an external field in the Ising model, say. 

Bosonizing the action gives for the partition function, 

Z^™™) = J V(t)VTTV^Vil:e-^ , (2.7) 

where 

S= [ (fx\ ^(x) (^ + M) i){x) + ^(t)'^{x) + ^-k'^{x) + (t){x)^[x)il;{x) + ■K{x)i){x)i-is'4^{x) \ , 

(2.8) 

where (pix) and 'ir{x) are auxiliary boson fields. The chiral transformations now represent rotations 
between these auxiliary fields. 



3 Lattice Regularization with Wilson Fermions 

Lattice regularization of the bosonized Gross-Neveu model, with Wilson fermions, leads to the 
action 

S'P^''[(/>,7r,'0,^/;] = 5j?^[^,V'] + 5'(int)[(A,7r,Vi,V'] + 5'(bosons)[^ , (3.1) 

where 



4 = ^ E Hn)Hn) - 2 E - 7M)V'(n + A) + + A)(l + 7/.)V'(n)} , (3.2) 

'S'(int) = E (t^{n)i^{n)i){n) + ^ T:{n)i){n)i-isi^{n) , (3.3) 

n n 

and 

5'(bosons) = ;7tE'?^^('^) + TtE^^("') • (3-'^) 

The lattice sites are labeled = —N/2, . . . , N/2 — 1, and is the number of sites in each of the 
two directions, which we assume to be even. Appropriate tuning of the two couplings and g"^ 
may allow recovery of chiral symmetry in the continuum limit (see [ pO| ] for discussions). 

Lattice Fourier transforms are defined as 

/H= (^)'E/»^''-"" ' /(p)=a2^/(n)e-^-- , (3.5) 

where 

= -^Pt^ ' (3-6) 

and where are integers or half integers depending on the field type and the boundary conditions. 
We henceforth drop the tilde on Fourier transformed field variables. 

The fermionic part of the action can be expressed in terms of momentum space variables as 



1 1 



a- N"^ 



Here M^^\q,p) are 2x2 matrices and 

M(^) {q, p) = M(0) (g, p) + M("^*) {q,p) , 



with 



M('\q,p) = 6p,,M^'Hp) , 



At At 



and 



M^'-'\q,p) = ^ E e*^''"'^"" [<^(^) + . 
Integration over the Grassmann variables gives the full partition function 

Z = J V(j)V7:VipV'il;exp (-5^7^) oc (detM^^^) oc (j[Xaip) 



(3.8) 

(3.9) 
(3.10) 

(3.11) 
(3.12) 



with Xa{p) the eigenvalues of the fermion matrix and the expectation values being taken over the 
bosonic fields. Note that there is no hopping parameter dependence in M^™*^. 



In the free fermion case the partition function is simply proportional to 

detM(o) = []AW(p) , 



where Xa^ (p) are the eigenvalue solutions of 



m(0)(p)|a(0)) = a(0)|A(°)) 



(3.13) 



(3.14) 



Using the representation ( p.2| ) for the Dirac 7-matrices, the solution to this problem is easily found 
to be 



|A(°)(p)) 



A(°)(P) 




1 

\a sinpia+i sinp2Q 



X cosp^a + i(-l)". 

Ai=l 



^sin^p^a 



(3.15) 



(3.16) 



where a = 1,2. These eigenf unctions form a complete orthonormal set. As is usual for Grassmann 
variables, we impose antiperiodic boundary conditions in the temporal (1-) direction and periodic 
boundary conditions in the spatial (2-) direction. With these mixed boundary conditions the 
momenta for the Fourier transformed fermion fields take the integer or half-integer values, 
p^ = -N/2 + 1/2, ... , N/2 - 1/2 and p2 = -N/2, N/2 - 1. Then, the eigenvalues (p^ in the 
free fermion case are either two-fold or four-fold degenerate, the former being the case if p2 = or 
-N/2. 

In the free fermion case the Lee- Yang Zeroes |23| are given by Xa\p) = 0. From ( 3.16| ), this is the 
case at 

1 

2k 



Vaip) = E cosp^a - i(-l) 
11=1 



. ^sinp^a 



(3.17) 



The lowest zeroes (with the smallest imaginary parts) correspond to 



p = {±{N/2-l/2),-N/2) , {±1/2, -N/2) , (±(A^/2 - 1/2), 0) and (±1/2,0) , (3.18) 

impacting onto the real axis at 1/2k = —2,0,0 and 2 respectively. These are precisely the three 
nadirs of the Aoki cusps in the Gross-Neveu model (see Fig. 1). 

Note that the zeroes in the upper half plane are given by a = 1, while their complex conjugates 
correspond to a = 2. Note, further, that the zeroes in (|3.17 ) are two or four fold degenerate in the 



momenta. I.e., these zeroes are invariant under —Pfi- This transformation is just a rotation 

through an angle vr in the space-time plane. The lowest zeroes ( ^.18 ), which are responsible for the 



critical behaviour of the free model, are actually two-fold degenerate. There is also a symmetry 
under pi ^ p2 which is manifest in the infinite volume limit. This is equivalent to a trivial rotation 
by 7r/2 in the pi-p2 plane, followed by reflection through the p2 axis. Since this reflection is through 
the spatial axis, this transformation is, in fact, parity. I.e., apart from rotations in space and time, 
the critical points are left unchanged under the parity transformation. 



4 A New Weak Coupling Expansion 

The usual weak coupling expansion of the full determinant for a general fermionic field theory is 
the Taylor expansion of 

detM(W) = detM(o) xdet(MW"'M(W)) 

= det M*^°^ exp tr In (i + m(°) ^m('°*)) . (4.1) 

This expansion is additive in nature and, from it, the ratio of full to free fermion determinants may 
be written, 

detM(W) ^ Mt^ 1 ^ Mjf'>M^'> l^Mth^^ 

detM(o) "^^^ a(°) 2.^^ a(0)^(0) +2.^^ ^^0)^(0) +••• • ^^.2) 

Here the indices i and j stand for the combination of Dirac index and momenta {a,p) which label 
fermionic matrix elements, so that M^^™*^ = Mj^^^^^^^^ represents {Xa\p)\M^^'^^\p, q)\X^^\q)) . The 
traces in (|4.2|) are, in fact, the diagrams which contribute to the vacuum polarization tensor. 



Setting 



ti = (Mt^) , (4.3) 
s^, = s,, = {Mt^Mff) , (4.4) 
= t,, = {Mjf^M^)-s,, , (4.5) 



the ratio of the interactive and free partition functions may be written 



detM(W)\ 2N^ , 2JV2 

^-l + ET|T-^ET7feuT + --- • (4-6) 



detM(o) ^ a(0) 2 X (0)^(0) 

1=1 i,j=i 



This Taylor expansion is analytic in 1/2k with poles at A^'^'' = or 1/2k = r]f\ 

The Wilson fermion matrix M^^^ is a 2N'^ dimensional square matrix given by ( |3.8| )- (|3.11| ). Its 
determinant, and the bosonic expectation value thereof, are therefore polynomials of degree 2N'^ 
with corresponding number of zeroes. As such, the latter may be written (up to an irrelevant 
constant) 

(detM(W)) = ]^(l/2^_^^) , (4.7) 

i=l 

where r/j represents rj^ip) and are the Lee- Yang zeroes of the interactive model. These are the 
quantities to be determined at weak coupling. 

Writing 

= r/i - r/i(0) , (4.8) 

gives, now, a new type of weak coupling expansion for the ratio of partitions functions, which is 
'multiplicative' rather than additive in form, 



detM(o) n( ^(0) ') n(i i/2K-,fj • ^'-'^ 



Note that the expression ( [4.9D , like its additive counterpart, (|4.6| ), analytic in 1/2k with poles at 

(0) 



Expanding (O) gives 



detM(W)\ . . ,2iV2 2JV2 . . 



Let {n} denote the n^^ degeneracy class in the free fermion case, so that the eigenvalues 
Xni = . . . = XnD„ are identical to Xn^ , say, with Dn = 2 or 4. Take the hopping parameter to be 
complex and arbitrarily close to a free fermion zero, 

l/2K = r/W+e . (4.11) 

The additive and multiplicative expressions (4.6) and ( [4.1C1| ) for the ratio of partition functions may 
now be expanded in e~^. Indeed, ( |4.6D gives 

detM(W)^ -2I -1/ tnj \ 



(4.12) 



while ( |4.10D yields 



detM(W)\ ^ ^ 

detM(o) ^ 2 



+.-'^ E E T55^^J+0(.») ^ (4,13) 



Equating these two expansions to 0(e ^) gives 



E 



An. A, 



riiTij 1 



(4.14) 



ni,nj(^\n\,ni^nj 



ni,rajg{n} 



while to 0{e-^) it gives 



E 1- E 



A, 



nie{n} \^ j^{n} V' 

having taken the e ^ hmit 
Let 



(0) __ JO) 



E 



JO) JO) 



E ' 



(4.15) 



U 
tij 
A, 



r?«+7?f) + 0(5^) , 



(4.16) 
(4.17) 
(4.18) 



where tj^'' and r/j^^ are the order g contributions to the expectation values of the matrix elements 
and zero shifts and where t^^\ t^-^^ and tjI"^^ are their order g^ equivalents. The 0{e~^) equation to 
order g is 

E — ~ E ) 



(4.19) 



niG{n} 



niG{n} 



and its order g'^ counterpart is 

E vi?^ 

ni£{n} 



ni€{n} 

Also, the 0{e~'^) equation, ( [4.14| ), is, now. 



.(2) , .{l).(l) 
'-rii ^ (0) (0) 



(4.20) 



E = E + E 

ni&{n\ ni,nj£{n\ \niG{n} 



With relations ( [4.19D -(4.21), the multiplicative expression ( 4.1C| ) recovers ( |4.6| ) to 0{g'^ 
equating the O (e*^) contributions to (4.12) and ( [4.13 ) yields no extra information. 



(4.21) 
Thus, 



The partition function zeroes are 'protocritical points' in the sense that they have the potential to 



become true critical points [24|. In the limit of infinite volume, the lowest zeroes impact on to the 
real hopping parameter axis precipitating the phase transition. The real parts of the lowest zeroes 
are therefore pseudocritical points in the statistical mechanics sense. 

In the free case, the lowest zeroes, and those responsible for criticality, are two fold degenerate. 
One expects critical behaviour in the weakly coupled case to be governed by their equivalents there. 
The two equations, ( [4.19| ) and ( 4.21D , allow full determination of the first order shifts to two- fold 
degenerate zeroes. Indeed, 



n 



(i) = l 



f(l) _ + Jft(l)+t(l)l2 , 4,(2) 



(4.22) 



where rii G {n} for i = 1 or 2. The second order equation, ( f4.20| ), in the two-fold degenerate case is 

/^(l) , ,{l)^^(l) ^(2) (2) 
„(2) , (2) _ _ ,(2) _ ,(2) , V- l^n, + In^ )l- + l-^^ + l-^^ 

'Ini ^ 'In2 — ''ni ''712 ^ (0) (0) " {'i.-^O) 

j<^{n} Vn -Vj 

To find the individual shifts, let 

^i? = Vi^^+S'^'^ , (4.24) 

^1? = ^i^)-'^^^^ • (4.25) 

J2) 



Their average, rjn , is determined directly from ( [4.23 ). Removing the expectation values over 



the bosonic fields converts the zeroes to the shifts in the eigenvalues of the fermion matrix in the 
presence of a small perturbation, M*^™*). The problem of determining such shifts is simply (two-fold 
degenerate) time independent perturbation theory. Indeed, one finds, for example, 

in which (^^^^^ = '^^^^^ This recovers time independent perturbation theory if 

1 A/f A/f A/f A/f 

^0 - 2 ^ x{0) ,(0) ' ^^-^^^ 

whence 

1 t(2) + 5(2) _ ^(2) _ ^(2) 

"2.^, (5^ • (^-2^) 

Finally, the full expression for the second order shift in an erstwhile two-fold degenerate zero is 

^^.(2) , .(2)^ , ^(1)^(1) , J2) J2) 

2 ~ 



j^{n} Vn Vj 



5 The Zeroes and Phase Diagram of the Gross-Neveu Model 



The interactive part of the fermion matrix (3.11) may be split into 



M(-*)(c7,p) = Mf*)(g,p) + M(-t)(5,p) , (5.1) 

where 

Mf''\q,p) = 4E^^^'"''^">('^)= (l^)''^(9-^') ' (5-2) 

M(-*)(g,p) = -l^e^(P-'?)-vr(n)i75= (^)'vr('Z-p)i75 . (5.3) 



One notes that the momentum dependency of the bosonic field variables involves even integers, 
so the bosons have periodic boundary conditions. The generic matrix elements required for the 



calculation of (O), (O) and (4.5) are 



. _^^E,,sinp^asing^a + z(sinpiosing2a - sinp2asingia) 



\NaJ ^'2 
^ ^ ^^^_^^E/,sinp^asing^a + z(sinpiasing2a - sinp2asingia) 

/J2f, sin^ p^,aJj2^, sin^ p^a 



In the (generalized) Gross-Neveu case, the pure bosonic action is given by (3.4). The pure bosonic 
expectation values in momentum space are thus 

im) = (vr(A:))=0 , (5.6) 

imH-k)) = N'^a^W^, , (5.7) 

{n{k)7T{-k)) = N''a^%l . (5.8) 

The bosonic expectation values of the matrix elements required in the calculation of the shifts 
(|42^ ) and ( p9| ) are then 



ti = t„,p = , (5.9) 
Sij = •S(c,,p)(^,5) = -j^ , (5.10) 

^ ^^-)(^^^)-^V^V-^) VE.sin^..E.sin^.. "7 • ^ ^ 

Prom these equations, together with ( [4.22| ) and ( f4.29| ), the 0{g) and 0{g^) shifts for the erstwhile 
two-fold degenerate zeroes, rya(±|pi|,P2) (for p2 = or —N/2), are, respectively. 



r?«(±|pi|,P2) = ±i , (5.12) 

So the two four-fermi interactions in fact contribute the same amounts to the shifts in the zeroes. 

In the thermodynamic limit these lowest zeroes become the true critical points of the theory and 
their determination amounts to determination of the phase diagram. I.e., the phase diagram is 
given to order g'^ by the limit 

-1^ = Jim U^\p) + 7?(^)(p) + r?(2)(p)} , (5.14) 



where p is the momentum corresponding to the lowest zeroes. The first order shift in ( 5.12| ) gives the 



relative separation of the erstwhile two-fold degenerate zeroes and vanishes in the infinite volume 
limit. The average shift is represented by (5.13) and is second order. The shift in the corresponding 
critical point is 



where 



hm 7]l^>{±\pi\,p2) = -{gl+gl) lim c{N) , (5.15) 

Af— >CXD W— >oo 



(/3,g)^{(Q,p)} Va (P) - ri/3 [Q) 

One finds, numerically, that the imaginary contribution to this factor vanishes in the thermody- 
namic limit, meaning that these zeroes indeed impact on to the real hopping parameter axis. The 



real part of (5.16) becomes an A^-independent constant whose actual value depends on the free 
zero from which it evolved. Indeed, ( ^.161 ) approaches approximately 0.77 and —0.77 for {pi,P2) = 
(±1/2,0) and (±(A^ — l)/2, —N/2) respectively. These correspond to the rightmost and leftmost 



critical lines (see Fig. 1). Also, ( 5.16|) is approximately 0.2 and —0.2 for {pi,P2) = {^{N — l)/2,0) 



and (±1/2, — A^/2) respectively. These give the two lines that generate the inner cusp. The situ- 
ation is summarized in Table |l] where the critical hopping parameters in the free and interacting 
cases and the momentum indices of the corresponding zeroes are listed. 

Table 1: The critical points and the momentum indices of their corresponding partition function 
zeroes in the free and weakly interacting cases. 



{Pl,P2) 


(±i(7V-l),-f) 


(±i-f) 


(±i(7V-l),0) 


(±io) 


1/2k,(0) 


-2 








2 


l/2Kc{g) 


-2 + 0.77(5^ + 5^) 


0-2(5| + 5^) 


-0.2(91 + gl) 


2-0.77(5^ + 5?) 



The actual phase diagram for weakly coupling is pictured in Fig. 1 for g^, = g-K = g (dark lines). 
The lighter curves are a schematic representation of the expected full phase diagram. One sees 
the degeneracy of the central free fermion critical point is lifted and two critical lines emerge in 
the presence of weak bosonic coupling. These are the lines corresponding to the central cusp in 
Aoki's phase diagram. From ( p. 15 ), the central cusp can only be made vanish at this order, in the 
unphysical situation of imaginary couplings. The Aoki phase does not yet emerge to 0{g'^) from 
the left- and rightmost free critical points. 

In the free case, the zeroes and hence the critical points are invariant under momentum inversion 
p^ — > —Pfj,, corresponding to a rotation in space-time. While this degeneracy is lifted at finite size 




1/(2 K) 

Figure 1: The phase diagram for the Gross-Neveu model in the weakly coupled region Qcf, = g-w = g 
(to 0{g'^)) (dark lines) and a schematic representation of the expected Aoki phase diagram (light 
curves). 



in the interacting case, it is recovered in the limit of infinite volume. In that limit, the free zeroes 
and critical points are also invariant under the parity transformation pi <-> P2- This is no longer 
the case in the presence of interactions. Indeed, the inner pair of critical lines are interchanged 
under parity. The overall phase structure, however, remains the same. 

The situation is similar to the two dimensional Potts model. There, the partition function is 
invariant under a duality transformation which exchanges the high and low temperature phases. 
The critical point is that which is invariant under that transformation. Here, the zeroes, and hence 
the partition function, are invariant under parity. The phase structure is also unchanged by parity. 
However, parity even and parity odd regions of the phase diagram are interchanged. 



6 Conclusions 



A new type of weak coupling expansion appropriate for Wilson fermionic lattice field theories has 
been developed. This expansion is multiplicative, but recovers the standard additive expansion. 
Its multiplicative form allows the Lee- Yang zeroes of the weakly coupled theory to be extracted in 
a natural way. These zeroes, are protocritical points, which, if they impact on to the real hopping 
parameter axis, precipitate a phase transition there. 

The expansion is applied to the single flavour lattice Gross-Neveu model to track the movement of 
zeroes and thereby the critical points in the presence of bosonic field variables. This model shares 
features with QCD, one of which is expected to be the existence of an Aoki phase. 

Using the new weak coupling expansion, a phase diagram is obtained in the weakly coupled region 
which is consistent with that of Aoki. The widths of the Aoki cusps are analytically determined to 
second order in the couplings. The central cusp cannot be tuned away for real physical couplings. 
The lateral cusps do not yet emerge at this order. This is the answer to the question posed by 



Creutz in |16] for the single flavour Gross-Neveu model. 

Finally, while the full phase structure is unaltered by a parity transformation, such an operation 
has the effect of exchanging the critical lines forming the inner Aoki cusp. 
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